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Abstract 

The paper represents the method for construction of the families of particular solu- 
^NJ ' tions to some new classes of (n + 1) dimensional nonlinear Partial Differential Equations 

(PDE). Method is based on the specific link between algebraic matrix equations and PDE. 
Admittable solutions depend on arbitrary functions of n variables. 

^ ; 1 Introduction 

o : 

Many different methods have been developed for analytical investigation of nonlinear PDE dur- 
ing last decades. Especially attractive are methods for study of so-called completely integrable 
systems. The particular interest to these equations is enhanced due to their wide range of ap- 
plication in physics. We emphasize different dressing methods, which are based on fundamental 
^ . properties of linear operators, either differential or integral: Zakharov-Shabat dressing method 
[1, 2], d-problem [3, 4, 5, 6], Sato theory [7, 8]. 

We suggest the method for construction of the families of particular solutions to some new 
classes of (n + 1) dimensional nonlinear PDE, n > 2. It is based on general properties of 
linear algebraic matrix equations. Essentially we develop some ideas represented in the ref . [7] 
and recently in the ref. [8, 9]. In general, for (n + 1) dimensional PDE this method supplies 
solutions, depending on the set of arbitrary functions of n variables. The represented method 
works also for classical (2 + l)-dimensional PDE integrable by the Inverse Scattering Technique 
(1ST). In this case our algorithm is similar to the algorithm represented in refs.[7, 10]. 

The structure of the paper is following. First, we discuss general algorithm relating lin- 
ear algebraic equation with nonlinear PDE. Then we show that these PDE are compatibility 
condition for some overdetermined linear system of equations having different structure in com- 
parison with linear system associated with completely integrable nonlinear PDE. We give an 
example of (2+l)-dimensional system which can not be described in frames of classical dressing 
methods. 
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2 General algorithm 



As mentioned above, our algorithm is based on the fundamental properties of linear matrix 
algebraic equation 

W = $, (1) 

where \I/ = {ipij} is N x N nondegenerate matrix, U and $ are N x M matrices M < N. 
Namely, the solution of this equation is unique, U = and consequently the homogeneous 

equation with matrix \& has only the trivial solution. Thus, if we find transformation T which 
maps the nonhomogeneous equation (1) into the homogeneous equation ^U{U) = 0, then 
U(U) = 0. 

Let us show that such transformations can be performed by means of differential operators 
having special structure. For this purpose let us introduce two types of additional parameters 
x = (xi, . . . , x n ) (n = dim(x)) and t — (ti, t 2 , . . . ) with the following systems: 

Vx^VBi + Qd, i = l,...,n (2) 

(Bi and C\ are constant N x N and M x N matrices respectively) and 

Mi* = 0, Mi$ = 0, Mi = d ti +L h (3) 

where Li are arbitrary linear differential operators having derivatives with respect to variables 
Xj and constant scalar coefficients, so that the system (2) is compatible with the system (3). 
For the sake of simplicity in this paper we use only one parameter t, omit subscripts in the 
eq.(3) and use n-dimensional Laplacian for L: 

n 

M = d t + Y, (4) 

fc=i 

Hereafter indexes i, j and k run values from 1 to n unless otherwise specified. 

Let us study compatibility conditions for the system (2) itself, which has the following form: 

(VBj + QC^Bi + $ x .Ci = {^B, + ^Ci))Bj + $ x .Cj. (5) 

Require that matrices Bi and C« satisfy two conditions: 

CjBi-CiBj = Q, BjBi-BiBj^O, i± j, (6) 

and matrices Cj have the following structure: Cj = [ Pj | 0m,n—r ] , R < M, where Pj are 
M x R matrices with rang R and 0a,b means Ax B zero matrix. Then equation (5) is reduced 
to the next one: 

$^•-$^ = 0, (7) 
which results in the first nonlinear equation for U owing to the eq.(l): 

(Bj + UCj)UPi + U Xj P t = (B t + XJCi)XJP 3 + U x .Pj. (8) 
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Let us show that another nonlinear matrix equation can be derived using operator M.. For 
this purpose we apply operator M. to both sides of the eq. (1) and use eqs. (2) and (3): 

n / n \ 

= M§ = (M*)U + VMU + 2 <Xk*x k U Xk = ^lMU + 2j2^(B k + UC k )U Xk = (9) 

k=l V k=l J 

Since det(^) ^ one has the second nonlinear equation for the matrix U: 

n 

U t + J2 a k (U XkXk + 2(B k + UC k )U Xk ) = 0. (10) 
k=i 

Having eqs. (6,8,10) one can derive the complete system of equations for elements of the matrix 
V composed of first R rows of the matrix U. First equation exists if n > 2. In this case we 
can write the matrix equation without operators Bi using any three equations (8) with pairs of 
indexes (J,k) and (k,i) and relations (6) (C/j = V Pi): 

T, P ^( U ^ k - U k Xj + [Uk,U j ]) =0, (11) 

perm 

where sum is over cycle permutation of indexes i, j and k. In general, this equation is not 
complete system for all elements of the matrix V. To complete it we derive additional matrix 
equation, which can be done for any n. Eliminate operators Bi from the eq. (8) (i — 1, j — 2) 
using eqs. (6) and (10): 

n n 

P i U jt -P j U it + 2 a k {P i U k -P k U i )U jxk -2 ^ a k (PjU k - P k Uj)U iXk + (12) 

n 

J2 ak { p Vi*u*u - p i u ***u + 2Pk{U ^ - - 2P k(Ui Xk Uj - U jxh Uij) = 

k=l 

In particular, if aj = 0, j > 1, a.\ = 1, then eq.(12) reduces in 

P 1 (U 2t + 2U lxiX2 - U 2xiXl + 2(C/ 2 C/ 1 ) :C1 " 2U lxi U 2 ) - (13) 
P 2 (U u + U lxiXl +2U 1 U lxi )^0. 

Thus equations (11,12,13) do not depend on both parameter iV (which characterizes dimensions 
of the matrices in the eq.(l)) and matrices Bi, i.e. N is arbitrary positive integer, Bi are 
arbitrary N x N matrices fitting relations (6). 

We point on two trivial reductions of the eq.(13). 

1. Let matrix A exist such that AP\ = and AP 2 = Ip> (In is R x R identity matrix). 
Multiplying eq.(13) by A from the right one receives matrix Burgers equation for U\. 

2. If R = M = 2 and 



Pi 



1 

1 



, v 



Ul Vi 

u 2 v 2 



(14) 



then eq.(13) is reduced to the next system: 

n - r Xl x 2 - 2rw XlX2 = 0, q t + q XlX2 + 2qw XlX2 = 0, w XlXl - w X2X2 = qr, (15) 
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where functions r and q are related with elements of the matrix V by the formulae 

U! = ^(r-q + 2w Xl ), u 2 = ^(r + q + 2w X2 ), (16) 
vi = -(-r - g + 2w X2 ), v 2 = -(-r + g + 2iu xi ). 

Eq.(15) becomes Devi-Stewartson equation (DS) after reduction r = ip, q = 4>, t\ = it, where 
i 2 = —1, bar means complex conjugated value. 

We will see in the next section, that the case R = M does always correspond to the classical 
completely integrable (2 + l)-dimensional systems. 

Arbitrary functions of variables Xj (j — 1, . . . , n) appear in the solution V due to the matrix 
function $, defined by the system (7). The number of arguments in the arbitrary functions 
as well as the number of these functions is defined by particular choice of the matrices Pj and 
dimension n of x-space. If n = 2, then one has at most iV arbitrary functions of two variables, 
see Example. In general, for n-dimensional x-space and R < M we are able to represent 
examples with N functions of n variables. If R = M then $ may depend at most on N x M 
arbitrary scalar functions of single variable, which is in accordance with [7]. Detailed discussion 
of this problem is left beyond the scope of this paper. 



2.1 On the operator representation of PDE 

We derive the overdetermined linear system of PDE with compatibility condition in the form 
of eqs.(ll) and (12). First, introduce arbitrary R x N matrix function R(A) of the additional 
parameter A. Multiply eqs.(8) and (10) by R(A) exp(r/) from the left and introduce function 
* = Re^U 

n / n \ 

V = Y,B k x k - \ Y^a k B 2 k \t. (17) 

k=l \k=l J 

We get after transformations: 

* x .Pi - * x .Pj = i>P l U J - VPjUi, (18) 

n 

i> t + Y,*k (** k x k + WPkV Xk ) = 0. (19) 
fc=l 

If R = M, t.e. all Pj are square nondegenerate matrices, then the system (18), (19) is equivalent 
to the classical M x M overdetermined linear system for correspondent (2+l)-dimensional 
integrable system. In fact, one can express all derivatives of ^ with respect to Xj, j > 1 
through the derivatives of with respect to x 1 using equation (18). Both equations (18) and 
(19) are M x M matrix equations for M x M matrix function Thus eq. (18) can be taken 
for the spectral problem while eq. (19) represents evolution part of the overdetermined linear 
system. For instance, if M — R — n — 2, ai — 1, a k — 0, k > 1, and Pi have the form (14), 
then the compatibility condition of the linear system (18,19) is given by the eqs.(15). 

In the case R < M situation is different. ^ is R x M matrix function, while (18) is R x R 
matrix equation. Thus it can not be taken for the spectral problem. Also, eq.(19) involves 
all derivatives which form operator M. So, we have (n + 1) dimensional equations. Below we 
represent example of (2 + l)-dimensional system of this type. 
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3 Example 



Let M = 3, R = 2, Q = 2, a 2 = a 3 = 0, cui = 1, TV = 3A; + 2, k = 1, 2, . . . . Everywhere indexes 
i and j take values 1 and 2. Let 
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, b j2 = diag(A,-, A,-, . . .), j = l, 2, 
fc-i 



1 -1 
-1 1 
-1 1 



A, 



1 -1 
1 -1 
-1 1 



$ = [0i 02 03], where 0^ are iV-dimensional columns. Then equation (7) can be written 
in the form 0i X2 = 2xi , (f) 2x2 = 3xi , i.e. 0i = S XlXl , 2 = S XlX2 , 3 = S X2X2 , where S 
is arbitrary function of variables x\ and a; 2 . In view of eq.(3), we can write for $: $ = 
f^c^ki, k 2 )[l k 2 /ki k 2 /kf] exp[kiXi+k 2 x 2 — kft] dkidk 2 , where c(ki, k 2 ) is arbitrary column of 
N elements. Function \& solves system (2,3) and can be represented in the form \& = (^ + ty)e ri , 
where 77 is given by the eq.(17), = d x ^C\e^ v and ^0 is arbitrary constant N x N matrix. 
Then eq.(l) can be solved for the function U, its first R rows define matrix V. Elements of this 
matrix satisfy the system (13) which gets the following form: 



uit + wixixi - 2u 2xiX2 + 2 v 2xiXl + 2mu lxi - 2(v 2 U!) Xl + 
2v 1 u 2xi - 2(w 2 u 2 ) xi + 2(u 2 ui) xi + 4:V 2 v 2xi 



Wit - Uxixi + 2uia. 1X2 + 2(lUiU2)xi + 2u!V lxi - 2u 2 Ui xl 

»it - ^lxixi + 2^1x1X2 - 2 wi«ixi + 4vi^i xi - 2w 2 7j la;i + 2(w 1 v 2 ) xi 
U2t + u 2xixi + 2u 2 u lxi + 2v 2 u 2xi = 0, v 2t + v 2xixi + 2u 2 v lxi + 2v 2 v 2xi 
w 2t - w 2xixi + 2u lxiX2 + 2v 2xiX2 - 2v lxixi + 2u 2 w lxi + 2v 2 w 2xi 











4 Conclusions 

The represented version of the dressing method serves for wide class of (n+l)-dimensional PDE. 
It supplies solutions depending on arbitrary functions of n variables provided R < M. In the 
last case solution depends on the functions of single variables. By construction, equations have 
infinite number of commuting flows and are compatibility conditions for some specific linear 
overdetermined systems, which is equivalent to the classical linear problem if only R = M. 

The work is supported by RFBR grants 01-01-00929 and 00-15-96007. Author thanks Prof. 
S.V.Manakov and Dr.Marikhin for useful discussions. 
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